Coherent propagation of waves in random media with weak nonhnearity 
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We develop a diagrammatic theory for transport of waves in disordered media with weak nonlin- 
earity. We first represent the solution of the nonlinear wave equation as a nonhnear Born series. 
From this, we construct nonlinear ladder and crossed diagrams for the average wave intensity. Then, 
we sum up the diagrammatic series completely, i.e. nonperturbatively in the strength of the non- 
linearity, and thereby obtain integral equations describing both nonlinear diffusive transport and 
coherent backscattering of the average intensity. As main result, we find that the nonhnearity sig- 
nificantly influences the magnitude of the coherent backscattering effect. Depending on the type of 
nonhnearity, coherent backscattering is either enhanced or suppressed, as compared to the linear 
case. 

PACS numbers: 04.30.Nk, 42.25.Dd, 42.65.-k 



I. INTRODUCTION 

Multiple scattering of waves in random media has been 
an active area of research for about hundred years. Of 
special interest is the role of interferences between differ- 
ent multiple scattering paths. For a long time, it was be- 
lieved that all interferences are destroyed by the disorder, 
such that the propagation of the average intensity follows 
a simple diffusion process. But in 1958, Anderson showed 
that interference effects can lead to a localization of wave- 
functions in certain random lattices, corresponding to a 
complete suppression of diffusion [I]. Later, it was recog- 
nized that even in weakly scattering media, interferences 
between amplitudes propagating along reversed scatter- 
ing paths may survive the ensemble average [2] , leading to 
the phenomena of weak localization (reduction of the dif- 
fusion constant) and coherent backscattering (enhance- 
ment of the average intensity scattered into the direction 
exactly opposite to the incident wave) [3, 4, 5j- 

Whereas these mesoscopic physics effects have been 
thoroughly studied for the case of linear wave equations, 
important, yet unresolved questions concern the impact 
of interactions or nonlinearities on the multiple scattering 
interferences. For example, does Anderson localization 
still persist in the presence of nonlinearities? Good can- 
didates for studying these questions appear to be Bose- 
Einstein condensates in disordered potentials [B]. In the 
mean field regime, the condensate is still described by a 
single coherent wave function following a nonlinear wave 
equation (Gross-Pitaevskii equation). Thus, the coher- 
ence of the condensate is preserved in presence of (not 
too strong) interactions. Similar nonlinear equations also 
describe propagation of light in disordered nonlinear me- 
dia [71 [S]. In contrast, the situation is quite different in 
the case of electronic transport j9j, where the interac- 
tions combined with finite temperature effects give rise 
to dephasing mechanisms, which in general destroy the 
disorder-induced coherent effects. 



As an important step in the theoretical description of 
coherent effects in nonlinear disordered systems, a the- 
ory for coherent backscattering in presence of nonhnear- 
ity has recently been developed [101 EI]- Most impor- 
tantly, it has been shown that the usual picture of co- 
herent backscattering resulting from interference between 
pairs of only two counterpropagating amplitudes breaks 
down in the nonlinear regime. As a consequence of non- 
linear coupling between different partial waves, nonlinear 
coherent backscattering must rather be interpreted as a 
multi-wave interference phenomenon |12l 113] . This fun- 
damental property was not taken into account in earlier 
works [3 US], which therefore concluded that nonhnear- 
ity would leave the enhancement factor in exact backscat- 
tering direction unchanged with respect to the linear case 
- in contradiction to our results presented below. 

The present paper provides a complete and detailed 
presentation of the diagrammatic theory sketched in jlOj , 
which fully takes into account the multi-wave interference 
character of nonlinear coherent backscattering. To per- 
form the disorder average, we will define nonlinear ladder 
and crossed diagrams, which describe diffusive transport 
and coherent backscattering of the average wave inten- 
sity in the dilute regime k£ ^ 1 (with wavenumber k 
and mean free path £). We are able to sum up the di- 
agrammatic series completely, i.e. nonperturbatively in 
the strength of the nonhnearity. As final result, we will 
see that even very weak nonlinearities do significantly 
affect the coherent backscattering enhancement factor, 
which, depending on the type of nonhnearity, may be 
either enhanced or reduced, as compared to the linear 
case. These predictions of our theory agree with recent 
numerical simulations [TOl [16] . 

The outline of the paper is as follows. In Sec.|TTj we de- 
scribe the general conditions under which the nonlinear 
diagrammatic theory can be applied. In short, these are 
the following: first, we consider a nonlinear equation for 
the complex wave amplitude E(r), which we assume to 
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possess a unique stationary and stable solution. In order 
to perform the disorder average, we require, second, the 
disordered medium to be dilute. Finally, we neglect non- 
linear effects which become relevant only at larger values 
of the nonlinearity. 

In the present paper, we will develop our theory on 
the basis of two different physical models, corresponding 
to homogeneous or discrete nonlinearity, respectively. In 
Sec. |III| we start by considering the discrete case, i.e. a 
collection of nonlinear point scatterers placed at random 
positions. We obtain the formal solution of the wave 
equation as a nonlinear Born series, and introduce di- 
agrams representing the individual terms of this series. 
Furthermore, we identify certain types of diagrams (lad- 
der and crossed diagrams) which survive the disorder av- 
erage over the random positions of scatterers. 

Among these, the ladder diagrams describe incoherent 
nonlinear transport of the average wave intensity. This is 
the subject of Sec. |TV] We analyze the statistical fluctu- 
ations of the wave intensity (which influence the average 
intensity, since the nonlinearity couples different statisti- 
cal moments), define refractive indices for scattered and 
incident waves, and finally arrive at a set of nonlinear in- 
tegral equations for the average intensity inside the scat- 
tering volume. 

In Sec.jV] we analyze the crossed diagrams, which de- 
scribe the coherent backscattering effect as interferential 
correction to incoherent transport. We identify the build- 
ing blocks for the nonlinear crossed diagrams, and point 
out the rules according to which these elements are con- 
nected to each other. This leads us to integral equations 
for the crossed intensity, which finally determine the co- 
herent backscattering enhancement factor. We will dis- 
cuss under which conditions the nonlinearity reduces or 
enhances the coherent backscattering effect, as compared 
to the linear case. 



In Sec. |VI| we consider a second physical model, 
namely a disordered collection of linear scatterers embed- 
ded into a homogeneous nonlinear Kerr medium. Equiv- 
alently, the same equation also describes Bose-Einstein 
condensates in disordered potentials (Gross-Pitaevskii 
equation). For this model, our theory can be applied 
in almost the same way as for the nonlinear scatterer 
model - in particular, the finally obtained average trans- 
port equations turn out to have the same form. We solve 
these equations numerically for different values of the 
Kerr nonlinearity constant a, and thereby confirm the 
general considerations of Sec. |V) 



Finally, Sec. VII concludes the paper. Some more 



technical calculations are relegated to two appendices. 
App. |X] contains the complete set of crossed transport 
equations, as they result from the exact summation of 
all crossed diagrams. In App.jS] we prove a general rela- 
tion between diffuson and cooperon cross section, which 
is useful in order to determine whether the nonlinearity 
decreases or enhances the coherent backscattering effect. 



II. GENERAL SETTING 

In this section, we will expose the general assumptions 
which we will need in Sees. [TTT||V| for developing the non- 
linear diagrammatic theory. 



A. Stationary state 

First, we assume that our system is described by a 
complex wave function E(y,t) e C, which fulfills a non- 
linear random wave equation. Furthermore, we assume 
that, for (almost) any realization of the disorder, the field 
E{r,t) approaches, for t oo, a stationary state of the 
form E{v.t) — e*'^*£'(r), where the frequency lu is deter- 
mined by the external, monochromatic source. Whereas 
for a linear system, this assumption implies no loss of gen- 
erality, since the time-dependent behavior for arbitrary 
(non-monochromatic) sources follows simply by Fourier 
transformation with respect to w, the same is not true 
in the nonlinear case. For example, it has been shown 
that solutions of the nonlinear Gross-Pitaevskii equation 
exhibit permanently time-dependent behavior even if the 
external source is purely monochromatic [17j . For a sys- 
tem of finite size, however, we expect the assumption of 
a monochromatic stationary solution to be valid if the 
nonlinearity is small enough, and if generation of higher 
harmonics can be neglected. Using the language of non- 
linear optics, the latter condition implies that all non- 
linear susceptibilities of even order (i.e. x^'^\ X^'^\ - • • ); 
which appear for noncentrosymmetric media |18j . van- 
ish. Apart from that, the polarizability of the nonlinear 
medium may be an arbitrary function of the local inten- 
sity |S(r)|2. 

Under these conditions, the field E{y) fulfills a station- 
ary, nonlinear wave equation, see for example Eq. ^ or 
Eq. (54) below. This wave equation will serve as the 



starting point for our diagrammatic theory, see Sec. |III[ 
Let us note however, that the diagrammatic theory in its 
present form is not able to predict the maximum strength 
of the nonlinearity up to which the stationary assump- 
tion remains valid. This must be checked by alterna- 
tive methods, for example direct numerical simulations. 
For the case of the homogeneous Kerr nonlinearity, see 
Eq. ( 54 1 , a theoretical prediction of the instability thresh- 



old was performed in [7] . Numerical simulations for non- 
linear point scatterers [12] and the 2D Gross-Pitaevskii 
equation [16] indicate that, in these cases, the instability 
threshold is large enough such that important nonlinear 
effects (i.e. modification or even inversion of the coherent 
backscattering effect) are observed in the stable regime, 
for which our theory applies. 



B. Dilute medium 

The second assumption required for our diagrammatic 
theory is the validity of a dilute medium approxima- 
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tion. When calculating the average scattered intensity, 
this assumption allows us to discard all diagrams ex- 
cept so-called ladder and crossed diagrams. From the 
ladder diagrams, we then obtain the diffuse background 
of the average scattered intensity, whereas the crossed 
diagrams yield the interference contribution leading to 
coherent backscattering. In the linear case, this dilute 
medium approximation is known to be valid if fc€ 3> 1, 
i.e. if the mean free path i is much larger than the wave- 
length A = 27r/fc divided by 27r. We now assume the 
phase randomization mechanism, which makes it possi- 
ble to neglect all diagrams except ladder and crossed di- 
agrams, to be equally effective also in the nonlinear case. 
In Sec. |IIIB| we will then define nonlinear generalizations 
of ladder and crossed diagrams, which we sum up in or- 



der to obtain nonlinear transport equations, see Sees. IV 
andlVl 

We note that, in principle, the dilute medium condition 
A;£ ^ 1 might be modified by the presence of the non- 
linearity. This must be checked for the specific model 
under consideration. For example, it can be shown that 
the diagrams describing scattering processes induced by 
the nonlinear index fluctuations a|i?(r)p in the case of 
the homogeneous Kerr nonlinearity, see Eq. ( 54 1 , are nei- 



ther of ladder nor crossed type. Since these processes de- 
pend on the strength a of the nonlinearity, the condition 
fc^ ^ 1 alone is not sufficient to guarantee their negligibil- 
ity. Instead, an additional assumption for a is required, 
see Sec. |VI A) On the other hand, for the nonlinear point 
scatterer model, see Sec. Ill A where the nonlinearity is 
concentrated inside the scatterers, we expect the dilute- 
ness condition to be fulfilled for arbitrary strength of the 
nonlinearity, if the average distance between the scatter- 
ers is large enough. 



non-perturbative in the strength a of the nonlinearity. 
As we will see later, even a very weak nonlinearity can 
have a large impact on the final result (i.e. the height 
of the coherent backscattering cone), if the size of the 
system (optical thickness) is large. In order to access the 
regime where perturbative expansions in powers of a do 
not rapidly converge (even if a is small), we will there- 
fore have to sum up our diagrammatic series completely, 
without restriction to certain orders of a. 



III. MODEL AND DIAGRAMS 

In Sec.|ll] we have discussed general conditions for our 
diagrammatic approach. In order to present the theory in 
a concrete form, we now specify the nonlinear wave equa- 
tion under consideration, and define a concrete model for 
the disorder. 

In this paper, we will consider the following two phys- 
ical models: a collection of nonlinear point scatterers, 
see Sec. 



Ill A 



and a collection of linear scatterers em- 
bedded in a homogeneous Kerr medium, see Sec. |VI A[ 
In the second model, the Kerr nonlinearity induces an 
intensity-dependent refractive index, which modifies the 
average wave propagation between two linear scattering 
events. Similarly, also in the first model, the nonlin- 
ear point scatterers form an effective nonlinear medium, 
but - in addition to that - also the scattering events are 
nonlinear. In this sense, the nonlinear point scatterer 
model is the more general one, and we will therefore use 
it as our starting point for the diagrammatic derivation 
of nonlinear transport equations. The modifications nec- 
essary for the homogeneous nonlinearity will be discussed 
in Sec. rvra 



C. Neglect of nonlinear self-action effects 

Finally, our last assumption concerns the neglect 
of nonlinear self-action effects (self-focusing, filamenta- 
tion,...) originating from the intensity-dependent refrac- 
tive index [18j. We assume those effects to be negligi- 
ble on the length scale of the mean free path i. More 
precisely, we approximate the optical path between two 
scattering events at ri and t^, with |ri — — ^, by 
a straight line, see Eq. (22 1, which also means that the 



spherical geometry of field radiated by a point source is 
preserved. This is justified if the intensity-dependent re- 
fractive index n(r) varies slowly on the length scale i.e. 
i?|Vn(r)| <C 1. In principle, these effects could be taken 
into account within the diagrammatic approach. The 
determination of the exact optical path required for this 
purpose, however, would enormously increase the com- 
putational effort. 

In summary, all the three above general conditions for 
the applicability of our diagrammatic approach are ful- 
filled in the case of dilute media and weak nonlinearity. 
Let us stress, however, that our approach is nevertheless 



A. Description of the model: nonlinear point 
scatterers 

We consider an assembly of N point-like scatterers lo- 
cated at randomly chosen positions r^, i = 1, . . . , iV inside 
a sample volume V illuminated by a plane wave k^. We 
assume the field radiated by each scatterer to be a nonlin- 
ear function f{Ei) of the local field Ei — E{ri). Since we 
neglect higher harmonics, the radiated field must attain 
the form 



fiE) = g{I)E 



(1) 



where / = EE* is the local intensity, and g{I) is propor- 
tional to the polarizability of the scatterers, which may 
be an arbitrary (smooth) function of /. This results in a 
set of nonlinear equations for the field at each scatterer: 



E, 



where 



gife|ri-rj| 

47r|r,; - r. l 



(2) 



(3) 
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with k = |ki,|, denotes the vacuum Green's function. 
Note that the incident plane wave is normalized to unit 
intensity, which can always be achieved by a suitable re- 
definition of g{I). 

In the case of conservative scatterers, the optical the- 
orem imposes the following condition: 



Im.g(/) = -|5(/)|^ 
which has a general solution of the form 



9(1) 



(4) 



(5) 



with phase shift 5{I) G K. Let us note, however, that 
our diagrammatic approach is valid for arbitrary polar- 
izability g{I). Hence, we will also consider amplifying 
or absorbing scatterers in the following, where Eqs. ( 4|5 ) 
are not fulfilled. 

Finally, the field Ejj measured by a far-field detector 
placed in direction ko (with \k.r> \ = fc) is given by 



3 — ikn-ri 

AttR 



(6) 



where R denotes the distance to the detector. From Ed, 
the intensity scattered into direction kc, measured as a 
dimensionless quantity (independent of R) , the so-called 
"bistatic coefficient" 1201, results as 



r(kD) 



47ri?2 
A 



\ED{\^Ti)? 



(7) 



where A denotes the transverse area (with respect to the 
incoming laser beam) of the sample volume V . 



B. Diagrammatic representation 

Our aim is to calculate the mean scattered intensity 
(r(kD)), averaged over the random positions of scatter- 
ers. For this purpose, we employ a method based on the 
diagrammatic representation of the iterative solution of 
the nonlinear multiple scattering equations 



1. A single nonlinear scattering event 

First, let us examine a single event, described by the 
function f{E) — g{I)E. Expanding the polarizability 
g{I) = 9nl^ in a Fourier series, we obtain the follow- 
ing sum over the order n of the nonlinearity: 



f{E)^Y.9nE''^\E*y 



(8) 



71=0 



The terms of this sum are represented by diagrams like 
Fig. [T| which shows, as an example, / for n = 2 and /* 



a) 



b) 




t 



Figure 1: Diagrammatic representation of noniinear scatter- 
ing events. Tfie feft half (a) represents a scattering event of 
nonfinear order n = 2, i.e. the term g2E''^{E*Y in the series 
Q for the scattered field f{E). The right half (b) repre- 
sents a compiex conjugate event of order n = 1, i.e. the term 
gi(E*)^E in the compiex conjugate series for f*{E). 



for n = 1. In Fig. [TJa), the outgoing solid arrow rep- 
resents the vacuum Green's function Gij pointing from 
position rj, where the scattering event / occurs, to po- 
sition Tj, where the next scattering event will occur in 
the subsequent iteration step, see below. The incoming 
solid and dashed arrows represent the arguments Ej and 
E* of the nonlinear function /, each of which may have 
undergone different previous scattering events. Due to 
the neglect of higher harmonics, the number of incoming 
solid arrows must exceed the number of dashed arrows 
by one. Of course, the opposite is true for the complex 
conjugate /*, Fig.jljb). 



2. Many scattering events: nonlinear Born series for the 
wave amplitude 

Using single scattering events as building blocks, we 
can now compose diagrams for the wave amplitude Ejj. 
Formally, a solution of Eq. ^ is obtained by iteration: 
we start with £'1"'' = 0, i = I, . . . N on the right hand 

side, and obtain E^^^ = exp(ikL -Ti) on the left hand side. 
This result again inserted on the right hand side yields 
and so on. In this way, the solution of Eq. ^ 
is obtained as a nonlinear Born series, i.e. a sum over 
sequences of single nonlinear scattering events. The indi- 
vidual terms of this sum are represented by diagrams of 
the form shown in Fig. [2] Here, all the arrows connecting 
two squares correspond to vacuum Green's functions Gij, 
whereas the remaining ones denote incident or outgoing 
waves. 

Please notice the following general structure of all non- 
linear field diagrams: Every full square (■) of nonlinear 
order n exhibits exactly one outgoing solid arrow, n in- 
coming dashed, and n -I- 1 incoming solid arrows. Like- 
wise, every open square (□) has one outgoing dashed 
arrow, n incoming solid, and n -I- 1 incoming dashed ar- 
rows. 

Note that, for simplicity, we have omitted in Fig. [2] 
any indices labelling the different scatterers. The reader 



Figure 2: Example of a nonlinear field diagram, obtained by 
iteration of the multiple scattering equation, Eq. (pi. The 
squares ■ (or □) represent scattering events / (or f). Solid 
(or dashed) arrows connecting two squares represent vacuum 
Green's functions dj (or G*_,), and arrows with bar the in- 
cident plane wave Icl- Finally, the half-circle symbolizes the 
detector placed in direction ko, see Eq. ||6|. 



Figure 3: Example of a nonlinear ladder diagram. It is ob- 
tained from two nonlinear field diagrams for Eo and E^, 
respectively (see Fig. [2|, by grouping together the nonlinear 
events ■ and □ such that every vacuum Green's function 
Gij (solid arrow) is accompanied by a complex conjugate G*j 



(dashed arrow). Scattering events (see Sec. IVAll occur at 
scatterer 1, 2, 5, and 7, where ■ and □ form a pair, whereas 
scatterer 3, 4, and 6 contribute to average propagation (see 
SeclrVBl). 



should simply imagine that every single event (I or □) 
occurs at some given scatterer i. It may be that two 
events occur at the same scatterer - only neighboring 
squares (which are connected by a single arrow) must 
correspond to different scatterers, since only couplings 
Gij between different scatterers i ^ j occur in Eq. (j2]). 
(Repeated scattering by the same scatterer is already in- 
cluded in the polarizability function g). 

We note that the diagrammatic representation of the 
nonlinear Born series has already been examined in a 
previous article [21] (for the special case of a x^"^^- 
nonlinearity) , which, however, does not make use of the 
dilute medium approximation, as we do in the following. 



3. Nonlinear ladder and crossed diagrams for the average 
wave intensity 

Next, we consider diagrams for the detected intensity 
r oc E£iE*jj, see Eq. ([t]), which consist of one diagram 
representing E^, and another one representing i?^, and 
perform the ensemble average over the random positions 
ri, . . . , Tat of the scatterers. 

For this purpose, we will employ the dilute medium 
approximation kH. ^ 1, where £ denotes the mean free 
path, see below. From the linear case, it is well known 
that - in leading order of l/(fc£) - only the so-called ladder 
and crossed diagrams contribute to the average inten- 
sity {E]jE*^). The reason why those diagrams survive 
the ensemble average is that, in the exact backscatter- 
ing direction, ItD = — kt, all the complex phase factors 
contained in the vacuum Green's functions Gij cancel 
each other. (In other directions ^ — kL, this is still 
true for the ladder, but not for the crossed diagrams.) 
This disappearance of phase factors is achieved by group- 
ing together solid and dashed arrows, such that for each 
vacuum Green's function Gij there exists a correspond- 
ing G*j or G*j which cancels the phase factor. Since 



Gij — Gji (reciprocity symmetry), the solid and dashed 
arrows may point in the same or in opposite direction, 
giving rise to ladder or crossed diagrams, respectively. 

In the following, we will now apply the same idea to 
our nonlinear diagrams. A ladder diagram is obtained 
if all the vacuum Green's functions occurring in the di- 
agram for EdE'^ are grouped together in pairs GijG*j 
pointing in the same direction. An example is shown in 
Fig. [3] Here, we adopt the convention that two squares 
(■ and □) drawn next to each other correspond to the 
same scatterer i (thus forming a scattering event at i, see 
below), whereas all other events occur at different scat- 
terers j i. Hence, the diagram shown in Fig.[3]involves 
7 different scatterers. Moreover, since the ladder dia- 
grams describe the average intensity, an integration over 
the positions ri of the scatterers is implicitly assumed. 
Thus, the example shown in Fig. [3] represents the follow- 
ing expression: 



G|2 1 ✓"Y* /- » /"f 

, I2I Lf23Lf25^43"^35<J65'-^76'-'75 



(9) 



where Af = N/V denotes the density of scatterers. Note 
that the phase factors exp(±ikL • ri,?) representing the 
waves incident on scatterers 1 and 7, and the waves 
exp(=F«kD -rs) emitted from scatterer 5 cancel each other, 
and are hence not included in Eq. ([9 1. 

The example shown in Fig. [3] is also useful to under- 
stand the difference between scattering events for the av- 
erage intensity, on the one hand, and events contributing 
to average field propagation, on the other hand. Scat- 
tering events always occur as a pair of two single events 
(■ and □), i.e., at scatterer 1,2,5, and 7 in Fig.js] where 
the wave amplitude E and its complex conjugate E* are 
scattered by the same scatterer, contributing thus to the 
propagation of the average intensity (EE*). In contrast, 
the events at 3, 4, and 6 influence the average amplitude. 
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a) 



^□1 




Figure 4: Example of a nonlinear crossed diagram. It is con- 
structed in the same way as the ladder diagrams, see Fig. [3] 
with the diflterence that pairs of arrows may also point in op- 
posite directions. The reversed pairs of arrows (thick lines) 
define the crossed scattering path. 



In linear diagrammatics, these events give rise to the av- 
erage Green's function (G(r,r')). However, since it is not 
evident how the concept of Green's functions can be gen- 
eralized to nonlinear equations, we will postpone a more 
detailed discussion to Sec. IIVBI 

In a similar way as the ladder diagrams, we can also 
construct crossed diagrams, see Fig.|4] Here, in contrast 
to the ladder diagrams, the two finally outgoing arrows 
originate from two different scatterers. Between those 
two scatterers, one can identify a 'crossed scattering path' 
by following the reversed arrow pairs. 

The mathematical expression for the diagram Fig. [4] 
reads: 



Gi2? G23G^5G:3G53G;6G;7G75e'^''^+''°' ""^ (10) 



Since Gij = Gji, the two expressions Eqs. ( 9|10 1 for the 
ladder and crossed diagram are identical in the exact 
backscattering direction Icd = — kL. This property is 
a consequence of the reciprocity symmetry. Unlike for 
the linear case, however, it does not imply that, in total, 
the crossed and ladder contributions to the backscatter- 
ing signal are identical. The reason is that, in general, 
a nonlinear ladder diagram has more than one crossed 
counterpart. Fig. |4] for example, just shows one of sev- 
eral different ways of reversing the scattering paths of 
Fig. [3] Indeed, reversing the direction of arrows in all 
possible ways (consistent with the general structure of 
nonlinear field diagrams, see Sec. IIIB2), we see that 



the outgoing solid arrow may originate either from scat- 
terer 1, 3, or 5, and the dashed one from 1, 5, or 7. This 
yields in total 2 different ladder and 7 different crossed 
diagrams, whose value (for ko — — kj^) is identical to 
Eq. ([9|. Hence, unlike to the linear case, nonlinear co- 
herent backscattering does not originate from interfer- 
ence between only two waves propagating along reversed 
scattering paths, but must be described as a multi-wave 
interference phenomenon. 

For this reason, the nonlinear crossed contribution can- 



□ I 



b) 




2 



A 



1 



Figure 5: Building blocks for nonlinear ladder diagrams. Di- 
agrams where ■ and □ occur as a pair (a ), con tribute to 
the average scattered intensity K, see Eqs. ( I5|16[ ), whereas 
(b) and (c) represent average propagation (G) and (G*), see 
Eq. pf. 



not simply be extracted from the ladder intensity by em- 
ploying the reciprocity symmetry, but must be calculated 
separately. In the following two sections, we will show 
how the ladder and crossed diagrammatic series can be 
summed up into a closed form, thereby obtaining nonlin- 
ear transport equations for the ladder and crossed inten- 
sities. 



IV. INCOHERENT TRANSPORT: SUMMATION 
OF LADDER DIAGRAMS 

In this section, we will express the sum of all ladder 
diagrams in terms of an integral equation for the aver- 
age intensity (/(r)). For this purpose, we first identify 
the building blocks out of which all nonlinear ladder dia- 
grams can be composed. These building blocks are shown 
in Fig. [5] As already mentioned in Sec. |IIIB"3l we distin- 
guish between scattering events for the average intensity 
represented by a pair of single events (I and □) as in 
Fig. |5|a), on the one hand, and average field propaga- 
tion. Fig. [5f^b,c), on the other one. 

As shown below, the average propagation events mod- 
ify the coherent part /c(r) of the intensity, which is de- 
fined as the square modulus of the average field ampli- 
tude, i.e. /c(r) = It represents the intensity 
originating directly from the incident field kL , and prop- 
agating until position r without being scattered into an- 
other direction k ^ kL . The remaining part of the aver- 
age intensity, called the diffuse intensity /(i(r) and defined 

by 



(/(r))=/,(r)-|-/,(r) 



(11) 



describes the part of the average intensity which has un- 
dergone at least one scattering event. 



A. Nonlinear scattering and diffuse intensity 

Thus, the diffuse intensity at point r inside the sample 
volume is obtained as the sum of the intensities radiated 
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by each of the individual scatterers at positions r': 

/,(r)= / dr'\{G{ry))fK{r') (12) 
Jv 

where K{r') describes the average intensity emitted from 
point r', and |(G(r,r'))p the subsequent propagation of 
the intensity from r' to r. 

1. Average scattered tntensity 

We first focus our attention on the average scattered 
intensity K(r'), which we now calculate by means of the 
nonlinear scattering diagrams, Fig. [5](a). 

The example shown in Fig. [sja) shows a pair of scatter- 
ing events / (■) with nonlinear order n ~ I, and /* (□) 
with nonlinear order m = 2 (see Sec. Ill B 1 1. To obtain 



the general result, we have to take the sum over n and m, 
respectively. Hence, the total number of incoming solid 
and dashed arrows in Fig.jsja) equals n + m+l. Accord- 
ing to the general recipe for constructing ladder diagrams 
(see Sec. IIIB3I, all arrows are grouped in pairs of solid 
and dashed arrows, respectively. Among these n + m+1 
arrow pairs, some represent the diffuse intensity Id, and 
the remaining ones the coherent intensity Ic- If we now 
label all incident arrows with a certain index, for example 
Eq for the coherent mode, and Ei, with i > 0, for a wave 
emitted by scatterer i, summing over all incident fields 
amounts to calculating the average value of / 



n+m+1 



as 



follows: 




n+m+1 



(13) 



- E ("^r^ j fc!^^^"^"' (14) 

with Ic — |i?oP and I^ — J2i>o 1-^*1^ (^^^ below). Within 
the ladder approximation, only such terms contribute in 
average to Eq. (13), where each amplitude Ei is accom- 
panied by its complex conjugate E* . Due to the large 
number of scatterers, we may furthermore neglect the 

^ 0, is 



cases where more than one arrow pair \Ei\ 
emitted by the same scatterer. This leads to Eq. (|14|, 
which can be interpreted as follows: among all n + m+l 
incoming arrow pairs, we select k solid and k dashed ar- 
rows contributing to Id, and take into account a factor 
kl for all possible pairings between them. 

Furthermore, note that we have identified in Eq. (14) 
the sum J2i>o l-^iP '^^ ^^'^ intensities radiated by each 
scatterer with the diffuse intensity Id given by Eq. (12). 



This is justified since, when forming a ladder diagram by 
connecting several building blocks like the one shown in 
Fig.jsJ^a) to each other, every outgoing arrow pair serves 
as an incoming arrow pair for the next building block. 



Finally, the average scattered intensity results as: 

C30 

Kild, Ic)=^ ff"5m (r+^+i) (15) 



with Af = N/V the density of scatterers, and 
given by Eq. (14). Since Id and Ic depend on r, also K 



is a function of r, as written in Eq. ( 12 1. Hence, Eq. ( 12 1 



forms a self-consistent equation for the diffuse intensity 
/rf(r). Solving this integral equation by iteration, we ex- 
plicitly obtain the sum of all ladder diagrams, i.e. all pos- 
sible connections of the building blocks shown in Fig. |5] 



2. Statistical properties of the wave intensity 

The statistical distribution of the wave intensity, which 
is expressed in terms of its statistical moments in 
Eq. ( 14 1 , can also be derived in an alternative way as fol- 
lows: we assume that real and imaginary part of the dif- 
fuse field Ed, defined hy Ed — E — (E), follow a Gaussian 



distribution with mean value zero and standard devia- 
tion -^Z Id/2. This so-called "chaotic" distribution for the 
diffuse field results from a random-phase approximation, 
i.e. by adding many partial waves (i.e. Ed = J2i>o ^i) 
with independent random phases, thereby neglecting any 
correlations between them [25]. In this sense, the calcula- 
tion of the average intensity by means of ladder diagrams 
is equivalent to neglecting correlations between different 
scatterers. This approximation is known to be well ful- 
filled in the case of a dilute medium {k£ ^ 1). 

The above considerations allow us to rewrite Eq. ( 15 1 
in a more compact form: 



K^J^{\f\^)^Af{l\g{I)\^) 



(16) 



where the average is taken with respect to the chaotic 
distribution for the diffuse field. Note that the result of 
this average depends only on Id and Ic = K^-)!^- For the 
explicit calculation of Eq. (16), it is useful to know the 



corresponding distribution function of the wave intensity 
/: 



P(/) 



-(/+/c)//<i 




(17) 



where Iq denotes the modified Bessel function of the first 
kind. In the absence of coherent intensity (i.e. deep 
inside the bulk), Eq. (17 1 reduces to the Rayleigh distri- 
bution P(/) 1/^=0 = cxp(—I/Id)/Id- On the other hand, 
for vanishing diffuse intensity (i.e., a scattering medium 
with extremely small optical thickness), the fluctuations 
vanish, such that P(/)|/^_^o ^ S{I — Ic)- 

Using the intensity distribution function, Eq. (17), 
the average scattered intensity, Eq. (16), results as: 
K{Id, Ic)=Af J dIP{I)I\g{I)\^. In the following, we will 
also need averages of the form {Eh{I)) , where h{I) may 
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be an arbitrary function of /. Here, the relation 

can be shown to result from the Gaussian diffuse field 
distribution. 



Nonlinear average propagation and coherent 
intensity 



In order to complete the transport equation (12), the 
coherent intensity Ic = | {E) ^ and the average Green's 
function (G'(ri, r2)) remain to be calculated. Both result 
from the diagram shown in Fig. [sjjb). 



1. Average Green's function and nonlinear refractive index 
for the diffuse intensity 

This diagram can be interpreted as follows: the scatter- 
ing event (■) modifies the propagation G12 represented 
by the solid arrows from 1 to 2. Integrating over the 
position ra of the scattering event, it turns out that, if 
kri2 3> 1 (dilute medium), only positions situated on the 
straight line connecting ri and r2 contribute to the inte- 
gral. (In other words: average propagation results from 
forward scattering.) The phase factors relevant for this 
stationary phase argument originate from the two vac- 
uum Green's functions G13 and G32, i.e. the two vertical 
solid arrows pointing from the starting point 1 and to the 
end point 2, respectively. In leading order of l/(fcri2), 
this integral yields: 



^f I dr3Gi3G32 



2k^ 



ikri2Gi2 



(19) 



The result on the right hand side should be compared 
with the expression ikri2{n — l)Gi2, which is obtained 
if one inserts a refractive index n in the exponent of the 
vacuum Green's function, Eq. (|3|, and expands the re- 
sulting expression in first order of n — 1. Hence, we can 
derive an expression for the refractive index n as follows: 
we sum over all possibilities to group the incoming ar- 
rows into pairs. Since there are m+1 incoming solid and 
TO dashed arrows, we obtain: 



For a conservative medium, i.e. for nonlinear scatterers 
of the form given by Eq. ^ , the extinction path length 
is equal to the mean free path. 



Since the intensity distribution function P{I), Eq. ( 17 1 , 
with respect to which the mean value (...) in Eq. (20 1 
is evaluated, depends on /c(r) and /d(r), the nonlinear 
refractive index attains a spatial dependence n(r). To 
obtain the average Green's function describing propaga- 
tion between two positions ri and r2, we now integrate 
the refractive index along the optical path (i.e. points of 
stationary phase) connecting ri and r2, which, as men- 
tioned in Sec. |II C"| we approximate by a straight line: 

exp (ik\ri - r2| dsn{rM 
<°'----»- ^ 4.|r.-r.| <^^> 

with — sri -f (1 — s)r2. By writing n in the exponent, 
we effectively sum up a series of diagrams with not only 
one, as in Fig. |5jb), but an arbitrary number of forward 
scattering events occurring between r and r'. 

From Eq. (22 1, we obtain |(G(ri,r2))P, which appears 
in the transport equation (12 1, as 



l(G(ri,r2))P = 



exp 



Ci-r2\J^ ds/£{rs) 



(47r|ri-r2|)2 



(23) 



The name "average Green's function" must be taken 
with some care, since, strictly speaking, nonlinear equa- 
tions like Eq. ^ do not have a well defined Green's func- 
tion. However, a Green's function can be defined for a 
linearized version of the wave equation ^ . The polariz- 
ability of a single scatterer in response to a small probe 
field AE is then proportional to (df/dE), which yields 
the refractive index, Eq. (20), and corresponding average 
Green's function, Eq. (22). The probe field interpreta- 



tion is also consistent with the observation that, due to 
the diluteness of the medium, the field emitted by any in- 
dividual scatterer is very weak compared to the incident 
field or the sum of the fields emitted by all other scatter- 
ers. In the following, we will hence continue speaking of 
"average Green's functions" when referring to the aver- 
age propagation between two scattering events. 



2. Refractive index for the coherent intensity 



1 = ^ £.9m(m+l)(/'") = ^ 



?TI = 



df 



2fc2 \dE 



with (/™) given by Eq. (14). The second equality results 
from taking the partial derivative of /, see Eq. ([s]) , with 
respect to E (where E and E* are treated as independent 
variables). The imaginary part of the refractive index 
defines the extinction path length I according to: 



1 



2fclm{n} 



The probe field argument, however, does not apply 

(20) 

for the coherent field {E{r)), which is not a weak field 
- at least close to the boundary of the medium, where it 
equals the incident field. In the diagrammatic calculation 
of Fig. ^h), this difference is reflected by the fact that 
the vertical solid arrow starting from point 1 in Fig. ^h) 
now originates from the coherent instead from a diffusive 
mode. This renders the combinatorial counting slightly 
different: we select I out of the m incoming dashed arrows 
and I out of the m-l- 1 solid arrows for the diffuse modes, 

(21) 

and group them in pairs in U possible ways. Thereby, 
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we obtain the refractive index for the coherent mode as 
follows: 



1 = 



m=0 1=0 

2fc2 (E) 



TO + 1 
I 



(24) 
(25) 



and 1/4 = 2Mm{nc}. Eq. (25 1 is derived from Eq. (24) 
with help of Eq. p^ . 

The fact that the refractive indices are different for 
diffuse and coherent waves should not be surprising: be- 
cause of the nonlinearity, the effective medium is modified 
by the propagating waves themselves, and, due to their 
different statistical properties, this is done in a different 
way by coherent and diffusive waves, respectively. From 
nonlinear optics in Kerr media, for example, it is known 
[T8] that a strong wave affects the refractive index of a 
weak wave of the same frequency twice as much as it 
affects its own refractive index, cf. Eqs. ( 58|59 ) below. 

Using the refractive index nc(r), the coherent field re- 
sults as: 



{E{r)) = exp [ ik J dz'nc{v,j) 



(26) 



where z denotes the distance from the surface of V to r, in 
the direction of the incident beam, and r^/ = r — z'ltL/fc. 
The corresponding coherent intensity reads: 

/,(r) - |(i?(r))|2 - exp 1^- j\z' / Uv,,)^ (27) 

Finally, let us remark that the above treatment of average 
propagation must be modified in the case where the in- 
cident wave is a superposition of several plane waves. In 
this case, additional coherent components are produced 
by four- wave mixing processes [18]. These are described 
by diagrams which are neither ladder nor crossed dia- 
grams, but nevertheless fulfill a phase-matching condi- 
tion, like for example exp[i(ki -I- k2) • r] = 1 in the case 
ki = — k2 of two opposite incident waves. 



a 



□ 



1 m 



□ 



\ it 




b) 



\ " if 



□ 



f) 



\ 



□ 

it 



Figure 6: Building blocks for crossed diagrams. By connecting 
these building blocks, we can construct any crossed scattering 
path, see the thick lines in Fig. H Diagrams (a- d) rep resent 



the crossed scattering events K^"-"', see Eqs. j33|36| , and 
(e,f) the crossed propagation events r, see Eq. (41 1, and r*, 
respectively. 



where Zq denotes the distance from r to the surface of 
V in the direction ko of the outgoing beam, and r^' = 
r + z'ko/k. 



COHERENT TRANSPORT: SUMMATION 
OF CROSSED DIAGRAMS 



The bistatic ladder coefficient, Eq. (28 1, expresses the 



incoherent part of the backscattered intensity, where all 
interferences are neglected. For a dilute medium, this 
is a good approximation if the detector is not placed in 
exact backscattering direction (jkn + kL| > l/£), but 
for ko — — kL, the coherent backscattering contribution 
resulting from the crossed diagrams must be taken into 
account. 



C. Incoherent component of the backscattered 
intensity 

We have now all ingredients at hand to determine 
the finally measured quantity, the intensity F(kD) scat- 
tered into direction kri, see Eq. ^. First, we solve the 
coupled system of integral equations, Eqs. ( 12|27| , in 
order to find the coherent and diffuse intensities /c(r) 
and /d(r). These, in turn, determine the intensity K(r) 
emitted from r, see Eq. (16 1, the refractive index n(r), 
see Eq. (20), and the extinction path length i{r), see 
Eq. (21). Finally, the bistatic ladder coefficient results 
as: 



rL(kD) 



dr 

AttA 



exp 



dzJ/£{r,,)]K{r) (28) 



A. Building blocks for crossed diagrams 

In analogy to the summation of ladder diagrams in 
Sec. |IV[ we first identify the individual building blocks 
out of which any crossed diagram is composed. These 
are shown in Fig. |6] For simplicity, we only display the 
reversed arrow pairs forming the crossed scattering path, 
which are drawn as thick lines in the example of Fig.[4j It 
is sufficient to concentrate on the crossed path, since the 
remaining elements, i.e., the thin lines in Fig. [4] all form 
ladder diagrams, and hence can be expressed in terms of 
the coherent and diffuse intensities /^^^(r), according to 
the results of the previous section. 
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1. Crossed scattering diagrams 

Like for the ladder building blocks, Fig. [5] we distin- 
guish between scattering and propagation diagrams. The 
crossed scattering diagrams, Fig. |6ja-d) are very similar 
to the ladder scattering diagram, Fig.[5]^a). The only dif- 
ference is that, in the ladder case, the two outgoing ar- 
rows are always grouped together, whereas in the crossed 
case, the outgoing arrows are grouped together with in- 
coming arrows. According to whether or not the incom- 
ing arrows are connected to the same scattering event 
as the corresponding outgoing arrow, we obtain the four 
possibilities displayed in Fig. [6|a-d). Remember that, 
as mentioned above, the crossed diagrams may be sup- 
plemented by an arbitrary number of incoming ladder 
pairs (thin lines in Fig. [4]). In particular, in the case of 
Fig. |4][b,c) and (d), at least one incoming ladder pair is 
required in order to obtain diagrams consistent with the 



Note that k^"'^ 



) , whereas and k^'^^ are real. 



general structure described in Sec. IIIB 2 according to 
which each nonlinear event exhibits one outgoing arrow 
and n -|- 1 or n incoming solid or dashed arrows. 

For the explicit calculation of diagrams Fig.[6](a-d), we 
first denote, as in the calculation of diagram Fig. |5ja), 
see Eq. ( |14| , the number of dashed arrows entering the 
filled square (■) by n (which determines the number of 
solid arrows entering this square as n-l- 1), and, likewise, 
the number of solid (or dashed) arrows entering the open 
square (□) by m (or m + 1). Then, one of the incoming 
solid arrows and one of the incoming dashed arrows must 
be chosen as partner for the outgoing dashed and solid 
arrow, respectively. According to which of the cases (a-d) 
we consider, this yields a factor n -f 1 or to for the solid 
arrows, and to -I- 1 or n for the dashed ones. After that, 
the number of remaining arrow pairs is decreased by one, 
i.e., 71 + TO instead of n + m-l- 1 in Eq. (151. Thereby, we 
obtain: 



rii-\-7i \ 



Ac) 



Ad) - 



rm+nx 



which, usine; Eq. (H, we rewrite as: 



Of df* 
dE dE* 

df df 
dE dE* 
df^dr_ 
dE dE* 
df* df 
dE dE* 



(29) 

(30) 

(31) 
(32) 

(33) 
(34) 
(35) 
(36) 



Let us remind again that the mean values (...) are taken 
with respect to the intensity distribution function dis- 
cussed in Sec. 



IV A 2 see Eq. ml 



As already discussed above in the case of average prop- 
agation, see Sec. IV B[ the expressions on the right hand 
side of Eqs. ( [33p6 r can be interpreted in terms of probe 
fields. Here, we need two independent probe fields AE 
and AE* , which represent the two reversed fields propa- 
gating along the crossed scattering path. The probe fields 
act as derivatives d/dE or d/dE* either on / or /*, what 
yields in total four different cases (a-d). Note that in the 
case of linear scatterers (/ = gE, with g = const.), cases 
(b-d) vanish, and only case (a) survives. 

To simplify the notation in the following equations, we 
define: 

,37) 

.^."> + .«^^(i^(ri|)) (38) 

We note that k — k* e M. 

If one of the incoming arrows participating in the 
crossed path originates from the coherent mode, the 
above calculation is modified in a similar way as for 
Eq. (24 1. As in Eq. (251, the result is that, in the case of 



a solid arrow, the derivative d/ (dE) is replaced by 1/ (E) : 



(E) 



dr 

dE* 



(E) 



df_ 
dE* 



(39) 



If the dashed arrow originates from the coherent mode, 
the complex conjugate expressions k* and k* must be 
used. The case where both solid and dashed arrows orig- 
inate from the coherent mode can be excluded, since the 
resulting diagram would be identical to a ladder diagram 
(single scattering). 



2. Crossed propagation diagrams 

Concerning nonlinear propagation, we can use the 
same average Green's function, Eq. (22 1, and refractive 
Eqs. ( [20p5t 



indices, Eqs. ( |20|25 l as in the ladder case, since aver- 
age propagation is an amplitude property, and symmet- 
ric against reversing the direction of propagation. To 
demonstrate this equality on the diagrammatic level, one 
just reverses the direction of the vertical "spectator ar- 
row" between 1 and 2 in diagrams Fig. [5]^b,c). However, 
it is also possible to reverse the direction of one of the 
arrows contributing to the "pump intensity" , i.e. the ar- 
rows coming from the left side in Fig. |5]^b,c). This yields 
two new diagrams, which are displayed in Fig.|6je,f). 

The physical interpretation of these diagrams is not 
straightforward: on one hand, they are very similar to the 
nonlinear propagation diagrams. Fig. [5]^b,c), from which 
they result just by exchanging two arrows. On the other 
hand, they can also be interpreted as nonlinear scattering 
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diagrams, since, as clearly evident from Fig. |6je,f), the 
crossed path changes its direction of propagation due to 
these events. Their mathematical evaluation, however, 
poses no difficulties. From Fig.[6][e), we see that we have 
to select two dashed and one solid arrow from all the 
arrows entering the nonlinear event. If the total number 
of incoming solid (or dashed) arrows is denoted by m (or 
m + 1), we obtain: 



a) 



■ h r oo 

2k ^ 

m— 1 



(40) 
(41) 



Diagram Fig. [6|(f) just yields the complex conjugate r*. 
The prefactor —i/{2k) in Eq. (41) results from the in- 
tegration over the position of the nonlinear event, see 
Eq. ( 19 1. If one of the incident fields originates from the 



coherent mode, d/dE (or d/dE*) is again replaced by 
1/{E) (or ll{E*)), i.e. 



2k 



(E) \{dE*)^ 



2k{E*) \dEdE* 



(42) 



Note also the following difference between the crossed 
scattering and propagation diagrams. Figs. |6ja-d) and 
(e,f), respectively: whereas in (a-d), only reversed arrow 
pairs are explicitly displayed, (e,f) also contain a ladder 
intensity entering the diagram from below (e) or above 
(f). In the exact version of the transport equations, see 
Eqs. ( A4| A8 ) , these will be represented by K{r'), i.e. 
the ladder intensity emitted from r'. In general, r' differs 
from the position r where the nonlinear event T(r) occurs, 
but if we assume that K{r') varies slowly on the length 
scale £, we may set r' = r. In this approximation, the 
contribution of the diagrams (e,f) to the cooperon cross 
sections will be given by T{r)K{r)£{r) and t* {r)K{r)£{r), 
see Eqs. ( |43|44| below. 



B. Connecting the building blocks 

1. Forbidden diagrams 

After having discussed the building blocks shown in 
Fig. [6] we now examine how they can be connected with 
each other in order to form a crossed scattering path. 
Like in the case of the ladder diagrams, see Sec. |IV| our 
aim is to derive an integral equation describing propaga- 
tion along the crossed scattering path through the scat- 
tering volume. We call the quantity which is transported 
along this path the "crossed intensity" C(r). To define 
the direction of propagation, we will choose the solid ar- 
rows in the following. Accordingly, we identify in each 
of the building blocks of Fig. [6] an "incoming" and "out- 
going" crossed intensity according to the direction of the 
solid arrow. 



b) 



A 



V 

A A 



it 



->n<r 

V 



Figure 7: Examples of forbidden (a) and allowed (b) combi- 
nations of the building blocks shown in Fig. [6] The forbidden 
combinations correspond to those forming a closed loop where 
two nonlinear events are connected by two Green's functions 
pointing in different directions. 



Similarly to the ladder case, see Sec.|IV] we could now 
obtain a crossed transport equation simply by identifying 
the incoming and outgoing crossed intensities, represent- 
ing both by the same function C(r). In this way, we 
would build crossed scattering paths by combining the 
individual building blocks of Fig. [6] in all possible ways: 
each outgoing crossed intensity may serve as an incoming 
crossed intensity for the next diagram. 

However, it turns out that some of the diagrams formed 
in this way exhibit closed loops. An example is shown 
in Fig. [Tj^a), where two scattering events (□ and ■) are 
connected by two Green's functions pointing in different 
directions. Such closed loops never occur in the series of 

Due to 
each 



2l E 



nonlinear field diagrams discussed in Sec. IIIB 2[ 
their construction as iterative solution of Eq. 
field incident on any given scattering event must form a 
valid nonlinear field diagram on its own. This property 
is obviously violated in Fig. |7][a), where none of the two 
scattering events can exist without the other one. Hence, 
we must avoid the occurrence of closed loops when con- 
necting the crossed building blocks with each other. 

As explained in the following, this is possible by split- 
ting the crossed intensity C(r) into two parts Ci(r) and 
C2(r). First, let us look at diagrams Fig. [6](a,c) and (e). 
Since these are the diagrams, where the outgoing crossed 
intensity is emitted by two different scattering events (I 
and □, respectively), we can connect them in all pos- 
sible ways without forming any closed scattering path. 
Thereby, we obtain the transport equation for Ci(r), see 
Eq. ( [46| in the following section. If one of the events (b,d) 
or (f) occurs, the crossed intensity changes from type Ci 
to C2, thereby indicating that the outgoing crossed in- 
tensity is now connected to a single scattering event (■). 
A closed scattering path would be formed if C2 was con- 
nected to a diagram where the incoming crossed intensity 
points towards the same event (□). This is the case for 
(c,d) and (e). Hence, the propagation of C2(r) may oc- 
cur only by the remaining diagrams (a,b) and (f). This 
yields the transport equation for C2(r), see Eq. ([47|. 
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2. Crossed transport equations 

As discussed above, the forbidden diagrams can be ex- 
cluded by splitting the crossed intensity into two parts 
Ci(r) and C2(r), respectively. The transport of Ci(r) is 
described by diagrams Fig. [6|^a,c) and (e), respectively. 
The quantities corresponding to diagrams (a) and (c ) ar e 
K^"^ and K,^'^\ whose sum is defined as k, see Eq. (37 1, 



whereas diagram (e) is approximately described by rKt, 
see Eq. (411 and the discussion after Eq. (42 1. Thereby, 



the three diagrams (a,c) and (e) governing the transport 
of Ci are summarized by the quantity 



CT(r) = ii{v) + t{v)K{v)1{v) 



(43) 



which we call "cooperon cross section" in the following. 
(More precisely, a corresponds to a cross section multi- 
plied by ^ttM , with M the density of scatterers. Thus, 
the dimension of a is one over length instead of length 
squared.) The transport of C2, performed by diagrams 
Fig.[6{a,b) and (f), is described by the complex conjugate 
cr*(r) of Eq. (III. 



Furthermore, the transition from Ci to C2 is deter- 
mined by diagrams Fig. [6]^b,d) and (f), or 



o-(r) = K(r) + T*(r)i^(r)^(r) 



(44) 



Note that ct = in the linear case. If the incoming or 
outgoing crossed intensity originates from the coherent 
mode, the above expressions are modified to CTc = Kc + 
TcKl and dc = kc + 'r*Ki. Finally, the "crossed coherent 
mode" Cc(r) is different from /c(r), Eq. ( pT] ), since only 
one of the two arrows (in our convention: the solid one) 
originates from the coherent mode (index ric), whereas 
the other one corresponds to an outgoing wave (index 
n). In total, we arrive at: 



Cc{v) 
Ci{r) 



expyikj^ dz'(nc(r^/) -n*(r^O)^ (45) 
dr'|(G(r,r'))P x 
(f7(r')Ci(r')+a,(r')C,(r')) (46) 
dr'|(G(r,r'))p(a*(r')C2(r')+ 
+a(r')Ci(r')+a,(r')Ce(r')) (47) 



from which the crossed bistatic coefficient (i.e. the height 
of the coherent backscattering cone) results as: 



Fr 



dr 

4:TtA 



X ((a;(r) + a:(r))Ci(r) + aUr)C2{r)) (48) 

The coherent backscattering enhancement factor is then 
defined as 77 = {Tc + i-e. the total backscat- 

tered intensity divided by the incoherent background. 



Writing the iterative solutions of Eqs. ( 46|47 1 formally 
as Ci = J2„{\G\^<Ty'\G\^<7cCe, and similarly for C2, it 
can be checked (using cr* — ct g K) that F^ G M, as it 
should be. 

In the expression for Cc(r), we have neglected events 
where the coherent mode is affected by a propagation 
event r^. Furthermore, Eqs. ( 45][48 ) are only valid in ex- 
act backscattering direction ko = — kL. An exact version 
of the crossed transport equations, where these assump- 
tions are relaxed, can be found in Appendix [A] 



Comparison between crossed and ladder 
intensity 



The interpretation of Eqs. (46 1 is as follows: we start 



with the crossed intensity Ci(?), or the coherent inten- 
sity Cc(r'), at point r' on the right hand side. Here, 
it undergoes a scattering event described by the corre- 
sponding cross section cr(r') or (Tc(r'), respectively. After 
being scattered, the crossed intensity propagates to r ac- 
cording to the average Green's function |(G(r,r'))p, see 



Eq. (22 1, which finally yields C'i(r) on the left hand side. 
A similar structure applies to Eq. ( 47 1 . 

For comparison with the ladder background F^, it is 
helpful to bring the ladder transport equation for /d(r), 
see Eq. (12l, into the same form as Eqs. ( 46|47 l. For this 



purpose, we write 



(49) 



with suitably defined "diffuson cross sections" a'^'^^ and 
ac^\ see Eq. (51) below. Then, Eq. (12 1 turns into the 
desired form 

/,(r) = / dr'|(G(r,r'))px 
Jv 

x(aW(r')/<i(r') + ^('')(r')/c(r')) (50) 

Making use of the intensity statistics defined by 
Eqs. p7p| , we can show that the following definition: 



(51) 



and, similarly, alf^ — Uc — d*^, fulfills the required prop- 



erty, Eq. (49 1. Eq. (511 can be rewritten as 



dl 



(52) 



In the conservative case, where the polarizabihty g{I) ful- 
fills the optical theorem, Eq. the diffuson cross sec- 
tions are related to the extinction path lengths through 
ct('*) = 4TT/i and a^/^ = 47r/4. 

Using Eq. (51), we can now discuss the relation be- 



tween the ladder and crossed intensities as follows: First, 
we note that in the linear case, where a = 0, diffuson 
and cooperon cross sections a^'^'^ and a are equal, as ex- 
pected from the reciprocity symmetry. This implies the 
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well known result that, in the linear case, the coherent 
backscattering cone height Tc equals the background T l 
(minus the single scattering contribution). 

Next, we turn to the case ct S K. Due to cr^'^'' S M, 
this implies also cr G K. As we can easily see from 
Eqs. ( |38|41|44| ), this case is realized for a purely imagi- 
nary nonlinear polarizability (g £ iM), or, in other words, 
if the nonlinearity modifies only the imaginary part o f th e 
refractive index. If a > 0, we find cr > cr('^) from Eq. m\, 



i.e., the cooperon cross section is larger than the diffuson 
cross section, and, consequently, Ci > Id- In addition, we 
have C2 > 0, and, consequently, also Ci + C2 > Id- We 
conclude that the nonlinearity enhances the cone height 
relative to the background if tr > 0. In a similar way, we 
can show that Tc/Tl is reduced if < 0. 

In order to find out whether a is positive or negative, 
we can use the following relation, which is valid in the 
absence of coherent intensity (i.e. deep inside the bulk): 



Re{ea}\j 



Idd{t(j^'^^) 
2" did 



(53) 



The proof of Eq. ( 53 1 is demonstrated in Appendix |Bj 
From Eq. ([53]) , we predict (in the case ct g M) that the co- 



herent backscattering cone is enhanced by the nonlinear- 
ity, if the scaled cross section increases as a function 
of the diffuse intensity, and vice versa. Note that the con- 
dition £o-('') — 47r just expresses the optical theorem valid 
for the case of a conservative medium. Thus, la'^'^'^ < Att 
corresponds to an absorbing, and la'^'^^ > An to an ampli- 
fying medium. Thereby, we arrive at the following con- 
clusion: A coherent backscattering enhancement factor 
larger than the linear value 2 (neglecting single scatter- 
ing) will be observed in the case of amplifying media, if 
the amplification increases with increasing wave intensity 
Id, or in the case of absorbing media, if the absorption de- 
creases with increasing Id (absorption saturation). In the 
opposite cases (gain saturation or increasing absorption) , 
the coherent backscattering enhancement factor will be 
diminished with respect to the linear case. This latter 
prediction of our theory was confirmed by exact numeri- 
cal simulations in [TU] . 

In the case of a general nonlinear polarizability g, 
which also (or only) affects the real part of the refrac- 
tive index, the cooperon cross sections a and a are com- 
plex numbers with non- vanishing phase. This implies 
the occurrence of phase differences between the reversed 
crossed paths, and, consequently, the height of the coher- 
ent backscattering cone is expected to decrease. 



VI. HOMOGENEOUS KERR NONLINEARITY 

The transport equations derived in Sees. |TV] and |V] 
are valid for the nonlinear point scatterer model with 
completely general form of the polarizability <?(/). In 
the present section, we will apply the diagrammatic 
method to a different type of nonlinear disordered sys- 



tem, namely linear scatterers embedded in a homoge- 
neous Kerr medium. In this case, the general structure of 
the transport equations will turn out to be the same as for 
the nonlinear point scatterer model - only the formulas 
for calculating the diffuson and cooperon cross sections 
are modified, as shown below. 



A. Model and transport equations 

We consider the following nonlinear Helmholtz equa- 
tion, describing propagation of light in a medium with 
randomly fluctuating refractive index ^e(r) and constant 
Kerr nonlinearity a: 

AE{r) + + 6e{r) + a\E{r)\'^'jE{r) = (54) 

We note that an equivalent equation (known as Gross- 
Pitaevskii equation) also describes propagation of Bose- 
Einstein condensates in disordered potentials |0. Con- 
cerning the disorder potential Se(r), we assume, for sim- 
plicity, short range correlations, i.e. {5e{r)6e{r')) ^ 
only if fc|r — r'| <C 1. Then, 6e{r) induces isotropic scat- 
tering, which is characterized by a single parameter, i.e. 
the mean free path Iq. 



The dilute medium approximation for Eq. ( 54 1 is valid 



if k£o ^ 1 and {al) Mq <C 1. The latter condition en- 
sures that scattering due to fluctuations of the nonlinear 
refractive index can be neglected 23] . This is automat- 
ically fulfilled if we assume that we are in the regime 
where Eq. (54 1 has a unique stable solution. According 
to |ff , this is the case (for a G M) if {alYb'^{Mo + 6) < 1, 
with b the optical thickness. 

To apply the diagrammatic method to Eq. (|54| 
first rewrite it as an integral equation: 



we 



E{v) = e'^^-^ + k^ J dr Go(r,r') x 
x(6e{r')+a\E{r')\^'^E{r') 



(55) 



which solves Eq. ( 54 1 for the case of an incoming plane 
wave. By iteration, we now obtain the nonlinear Born 
series in a similar way as for Eq. ([2]). When performing 
the ensemble average, some complications arise concern- 
ing the nonlinear scatteri ng e vents represented by the 
term a\E{r)\'^E{r) in Eq. (55 1. Since their position r is 



not restricted to a discrete set of dilute point scatter- 
ers, we must, in principle, take into account correlations 
of the form \a\^\{\E{ri)\'^E{ri)\E{r2)\^E*{r2)) with n 
close to r2. It turns out, however, that we may neglect 
these terms if we assume (a/)^Ho li as mentioned 
above. Since only linear scattering events remain, we im- 
mediately get the following results valid for linear point 
scatterers: 



rid) 



47r 




(56) 
(57) 
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Concerning nonlinear average propagation, we recover 
the same equations ( 20|24 ) for the refractive indices and 
Eq. (41 1 for the crossed propagation diagram as for the 
nonlinear point scatterer, if we replace, in these equa- 
tions, the density of scatterers times the polarizability 
by the corresponding expression for the Kerr nonlin- 
earity, i.e. Afg{I) —> k^al. Adding the linear index 
no — 1 = i/(2k£o), we obtain in total: 



n-l = + + 



1 = 



2Hn 



(58) 
(59) 



with corresponding extinction path lengths £ — 
l/(2fclmri) and £c — l/(2fclmric), and 

T = Tc = —ika* (60) 



The above expressions, Eqs. ( 56p0 1, can now be inserted 
into the transport equations for the ladder and crossed in- 
tensity, Eqs. ( 12|27|) a nd Eqs. ( [A4]|A8| . Also the approx- 
imate form, Eqs. ( 45]|47 1, as well as the relation (53) be- 
tween diffuson and cooperon cross section remain equally 
valid for the homogeneous nonlinear medium, with the 
cooperon cross sections a and a given by Eqs. ( 43|44 ). 



B. Results for conservative nonlinearity {a G R) 

First, we examine the case a G M, where energy 
is conserved (more precisely: the intensity flux j = 
Ke{iEyE*} fulfills the conservation law Vj — 0). In 
this case, we see from Eq. (56 1 and Eqs. ( 58|59[ ) that the 
quantities K (or, equivalently ct*^**^ and oi^^ £ and 



which enter into the ladder transport equations ( 12|27 1, 
are unaffected by the nonlinearity a. Hence, the average 
diffuse and coherent intensities /(i(r) and /c(r), as well as 
the background component of the scattered intensity, 
see Eq. (28 1, are the same as in the linear case. This is 
not surprising, since for a € M, the nonlinearity modi- 
fies only the real part of the refractive index, which, in 
the case of a dilute medium, has a negligible influence on 
incoherent wave transport. 

However, the nonlinearity does affect interference 
properties like coherent backscattering, which are de- 
scribed by crossed diagrams. To see this, we first deduce 
from Eqs. ( 56]|60 1 the cooperon cross section for a £ M, 
and obtain: 



cr(r) = crc(r) 
5-(r) = CTc(r) 



47r 

4 



[1 - zHoa (4(r) + /d(r))] , (61) 



inika {Ic{r) + Id{r)) 



(62) 



Of special interest is the fact that Rea ~ 0. From 
Eq. (53 1, we see that this is generally true (at least in- 
side the bulk where Ic = 0) for conservative media, where 
fcr'''' = An. Using the condition Retr = 0, it can be shown 



(again with help of the iterative solutions for Ci and C2) 
that: 



c 



Re 



dr 



Ct{v)a{v)C,{v) 



(63) 



tion 
Eq. 



If Cc e M (i.e. if the difference between n(r) and ric(r) is 
neglected) , Eq. ( 63 1 can also be shown in a simpler way, 
since then ReC2 = and ImC2 = — ImCi follow directly 
by splitting Eqs. ( 46|47 ) into real and imaginary parts. 

Due to Eq. ( 63 1 , solving the transport equation for 
Ci, Eq. (46), is sufficient for calculating the coherent 
backscattering enhancement. As already mentioned, this 
equation is very similar to the ladder transport equa- 
Eq. (50), with the difference that the cross section, 
(56), is a complex quantity a — \a\ exp^iip) . This 
can be interpreted as follows: at each scattering event, 
the nonlinearity introduces a certain phase difference 4> 
between the two reversed paths, which depends on the av- 
erage intensity (/(r)) = /^(r) -f /c(r) multiplied by k£Qa. 
As a result of this dephasing mechanism, the height of 
the coherent backscattering cone decreases with respect 
to the linear case. Note that the sign of a plays no role, 
i.e. we obtain the same value of Tc for focussing (a > 0) 
and defocussing (a < 0) nonlinearity, respectively. This 
also results in an initial decrease quadratic in a. 

In the conservative case, we hence recover the famil- 
iar picture of coherent backscattering as an interference 
between only two reversed amplitudes - however with a 
phase difference induced by the nonlinearity. At flrst 
sight, it might be surprising that the nonlinearity does 
not affect both counterpropagating amplitudes in the 
same way. This can be understood only by taking into 
account the many-wave interference character of nonlin- 
ear coherent backscattering. In the conservative case, 
every nonlinear event gives rise to a purely imaginary 
factor, i.e. t = —ika or r* = -|-ifca, such that most of 
the many interference terms cancel each other. At the 
end - as a consequence of the forbidden diagrams, see 
Sec. VBl - it turns out that only a single interference 
term remains, where the nonlinearity modifies only one 
of the two counterpropagating amplitudes "0 or ijj* , but 
not the other one. 

In the limit akiol ^ 1, we find from Eq. ( [M] ) the 
phase difference A(/) = Mk£oaI for a path undergoing M 
scattering events. In the case of a slab geometry (finite 
length L — Mq in z direction, infinite extension in x and 
y direction) the mean number (M) of scattering events 
for diffuse reflection is proportional to the optical thick- 
ness b of the scattering medium (in contrast to diffuse 
transmission, where (M) oc b^). From this, we predict a 
significant reduction of the coherent backscattering peak 
if bk£oaI ~ 1. Note that, if b and k£Q are large, already 
a very tiny nonlinearity is sufficient to observe a large 
impact on the CBS peak, which, in particular, still lies 
in the regime {aI)^b'^{k£o + &) < 1, where, according to 
[7], the wave equation has a unique stable solution. 

The above conclusions are confirmed by Fig. [8] show- 
ing the results obtained from the exact transport equa- 
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Figure 8: Coherent backscattering cone height Tc as a func- 
tion of the nonhnearity strength kloa for a slab with opti- 
cal thickness 6 = 2 (solid line) and 6 = 10 (dashed line). 
The initial decrease of the cone height due to dephasing be- 
tween reversed paths is faster for larger optical thickness. At 
large nonlinearities, the cone height Tc assumes negative val- 
ues, indicating destructive interference. The thick lines result 
from the complete summation of crossed diagrams, see Ap- 
pendix [X] whereas the thin dotted lines correspond to the 
simplified version, Eqs. ( |45|48[ ). 



tions, Eqs. (A4A9I, compared to their approximate 
form, Eqs. (45] 48 1, for a slab geometry with optical thick- 
ness 6 = 2 and 6 = 10, respectively. As expected, the ap- 
proximate form agrees very well with the exact result for 
large optical thickness and small nonlinearity. Further- 
more, the initial decrease, quadratic in a, of the cone 
height occurs much faster for 6 = 10 than for 6 = 2, as 
predicted above. Remarkably, the cone height assumes, 
in both cases, negative values for large nonlinearities, cor- 
responding to destructive interference in backscattering 
direction. This property has recently been observed also 
in numerical studies of the 2D Gross-Pit aevskii equation 



C. Results for absorbing/amplifying nonlinearity 

(a G m.) 

Next, we discuss the case (3 = ikloa, with /3 G M. 
Then, p > corresponds to an amplifying, and f3 < 
to an absorbing medium. According to Eqs. ( 58|59 1, we 
find for the extinction path lengths: 



1 



1 

1 



(1- 2/3(4 + /rf)) 
(l-/3(/, + 2/rf)) 



(64) 
(65) 



Note that, for /3 > sufficiently large, the extinction 
path length may, in principle, become negative (as a con- 
sequence of strong amplification). However, as we will 
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Figure 9: Background Fl (dashed lines) and cone height Fc 
(exact version: solid lines, simplified version: dotted lines) as 
a function of /3 = ikloa (amplifying nonlinearity), for a slab 
with optical thickness a) b = 2 and b) 6 = 10. The cone height 
is more strongly affected by nonlinear amplification than the 
background. Close to the critical value /3cr of the nonlinearity, 
the cone height exceeds the background intensity by more 
than a factor 3. 



find below, stationary solutions of the nonlinear trans- 
port equations typically exist only for small f3. Therefore, 
we assume /? <C 1, and hence ^ > in the following. 

By inserting Eqs. ( 56|57|60" l into the definitions of a 
and a, Eqs. ( 43|44 1, we obtain: 



cr = (Jr. = 



An 1 - 13(1, + h) 
£o 1 - 2/?(/e + Id) 
4tt I3{h + Id) 
lo 1 - 2/3(4 + Id) 



(66) 
(67) 



From the sign of tr, we can now conclude 
discussion in Sec. |VB 2 
creases for /? > 0, and decreases 



following the 
that the CBS cone height in- 
for /3 < 0. This is 



confirmed by Fig. |9] where the background intensity F/, 
and cone height Fc are plotted as a function of for a 
slab with optical thickness b — 2 and 6 = 10, respectively. 
Again, the exact solutions for Fp (solid lines) are com- 
pared with the simplified version (dotted lines), which 
fits well in the case of large optical thickness and small 
nonlinearity, as expected. More importantly, however, 
we see that the cone height is more strongly amplified by 
the nonlinearity than the background, such that F^ may 
exceed Fl by more than a factor three. This is a direct 
evidence for the multi-wave interference character of CBS 
in the nonlinear regime, as mentioned in Sec. |IIIB 3| 

However, we note that the backscattering cone height 
Tc does not become arbitrarily large. When increas- 
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ing the amplifying nonlinearity beyond a certain crit- 
ical value Pcv, we find that, for f3 > Pen the station- 
ary nonlinear ladder transport equations ( 12|27 1 have no 
solution. This behavior can be understood if we ap- 
proximate Eq. ( |12[ ) by a nonlinear diffusion equation. 
For this purpose, we assume that the average intensity 
I{z) = Ic{z) + Id{z) varies slowly on the length scale of 
a mean free path, i.e. IqI'{z) <^ 1 and £q/"(z) ^ 1. Ad- 
ditionally, we assume /3 ^ 1. Expanding Eq. (12 1 up to 



lowest non-vanishing order in those small parameters, we 
obtain: 



I" + 2(31^ 



-zjla 







(68) 




where we have approximated Ic{z) ~ exp(—z/£o). The 
boundary conditions read: /(O) — zo^'(O) = and I{L) + 
zqF' {L) — 0, with L = Mo the thickness of the slab, 
and the extrapolation length zq ~ 0.71i?o PT- From the 
solution of Eq. ( 68 ) , the background intensity is obtained 



Figure 10: Background intensity as a function of /? = 
ikioa (amplifying nonlinearity), obtained from the diffusion 
equation (681 for a slab with optical thickness b = 10. For 



as Fi = /p dzI{z)eyi\){~z/io). 



For small /? > 0, the nonlinear diffusion equation (68) 



has two different solutions, see Fig. 10 One of them (The 



solid line) reproduces approximately the iterative solu- 
tion of the ladder transport equations shown in Fig. |9] 
On the other hand, the iteration of the ladder transport 
equations diverges if the initial condition is chosen close 
to the other solution (dashed line in Fig. [lO| ). We sus- 
pect that the lower branch (solid line) corresponds to a 
dynamically stable and the upper branch (dashed line) to 
a dynamically unstable solution of the nonlinear diffusion 
equation. 

Finally, for (3 > /3cr, the nonlinear term 2/3/^ in the 
diffusion equation ( 68 ) enforces a too large negative cur- 
vature, I"{z) < 0, such that no solution exists which is 
positive, I{z) > 0, in the entire interval < z < L, 
and which satisfies the boundary conditions mentioned 
above. Hence, in this regime, our assumption that the 
wave equation ( 54 1 possesses a stable stationary solution 
(see Sec. II A I cannot be fulfilled. Instead, an explic- 



itly time-dependent treatment of the scattering problem 
would be required. 



VII. CONCLUSION 

We have presented a diagrammatic theory of nonlinear 
transport and coherent backscattering for dilute, disor- 
dered media with weak nonlinearity. The theory was ap- 
plied to two different models: nonlinear point scatterers 
and linear scatterers embedded in a homogeneous non- 
linear medium. The relevant nonlinear effects captured 
by the theory are the following: first, scattering cross 
section, refractive index, and mean free path depend on 
the average wave intensity and its statistical fluctuations. 
Furthermore, the refractive indices are different for coher- 
ent and diffuse waves, respectively. These properties are 
relevant for incoherent transport of the average wave in- 
tensity. Concerning properties of coherent transport, the 



< /3 < /3cr, two solutions are found. The lower branch (solid 
line) reproduces the iterative solution of the ladder transport 
equations ( |12|27| l, cf. the dashed line in Fig. [9|b), whereas 
the upper branch (dashed line) corresponds to a physically 
unstable solution. No positive solution of Eq. (681 exists for 



most important conclusion is that, due to nonlinear cou- 
pling between different partial waves, coherent backscat- 
tering turns into a multi-wave interference effect - instead 
of two-wave interference, as in the linear regime. 

How the multi-wave interference character affects co- 
herent backscattering quantitatively, depends, of course, 
on the phases between the various interfering amplitudes. 
Complete constructive interference - and hence nonlin- 
ear enhancement of the coherent backscattering effect 
- is found for the case of amplifying nonlinearity. For 
a conservative nonlinearity, we find effectively only two 
interfering amplitudes, however with a phase difference 
proportional to the nonlinearity strength times the aver- 
age intensity integrated along the whole scattering path. 
In this case, the coherent backscattering is suppressed (or 
even inverted) by the nonlinearity. 

The nonlinear enhancement of coherent backscattering 
predicted in the case of amplifying nonlinearity opens in- 
teresting perspectives concerning the phenomenon of An- 
derson localization, which (in the linear case) is expected 
to occur in the regime fc£ < 1 of strongly scattering me- 
dia. Is Anderson localization affected in a similar way 
as coherent backscattering, i.e., can it be suppressed or 
enhanced, depending on the type of nonlinearity? This 
is an interesting problem to be investigated in future re- 
search. 
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Cc{r) = 



Ci(r) 



K: 



■>□<■ 

Ci, C'c 



Figure 11: Graphical representation of the transport equa- 
tions for Ci and Cc, Eqs. I A4|A5| |. The letters 'a' and 'c' 
refer to the crossed scattering building block from Fig.|6] i.e. 
they represent the quantity Eq. (|37|. 



the DFG. 



Appendix A: EXACT SUMMATION OF CROSSED 
DIAGRAMS 



In the approximate version of the crossed transport 
equations, Eqs. ( 45p7 l, the crossed scattering and prop- 
agation diagrams, Figs. [6|a-d) and (e,f), are treated on 
an equal footing, such that their combined effect is sim- 
ply expressed by the cooperon cross sections cr and a, 
Eqs. ( 43|44 l. In this appendix, we demonstrate the ex- 
act calculation of the nonlinear propagation diagrams, 
and examine the conditions under which the approximate 
version is expected to be valid. 



1. Crossed propagation diagram 

Let us first look at diagram Fig. [6|e), which corre- 
sponds to the quantity t, Eq. (41 1. Similarly as in the 
ladder 



case. Fig. [5|^b), we denote the positions corre- 
sponding to the bottom and the top part of the diagram, 
which are connected by the vertical solid arrow, as ri 
and r2 , respectively. The position of the scattering event 
□ is denoted as Vs- Due to the condition of stationary 
phase, see Eq. (19 1, rs must be situated on the straight 
line between ri and r2, i.e. rs — sri + (1 — s)r2, with 
< s < 1. The pair of reversed arrows on the right hand 



side represents an incoming crossed intensity, which we 
denote by Cin(r3). From ri (the bottom), a ladder in- 
tensity is emitted, which is described by K{ri). In total, 
the diagram Fig. [6|e) represents the following expression 
for the crossed intensity at 

C(^)(r2) = / drii^(ri)|(Gi2)P x 
Jv 

xri2 / dsr(r3)Cin(r3) = (Al) 
Jo 

= [ dr3T(r3)an(r3)|(G23)P X 

Jv 

X / dpie-P'/'^Kin) (A2) 
Jq 

where ri = rs -I- pi{rs — r2)/r32, and the upper 
limit pj™^''-' is determined by the condition that ri 



must lie inside the sample volume V. Eq. (A2| re 



suits from Eq. (Al) by a change of variables (ri,s) 



(r3,pi). In the last line, we use a short-hand notation 

exp(-pi/^) := exp (^-pi dt/i^ra - trs + tri)j for the 

average of l/^(r) over the line of length pi connecting ri 
and r3. 



2. Equations for Ci 



According to our discussion in Sec lVBTj we recall that 
diagram Fig. ^e) contributes to the crossed intensity of 
type C'l . The remaining diagrams contributing to Ci are 
Fig.|6ja) and (c), which are described by the quantity k, 
see Eq. (p7|), as follows: 



C 



(a+c) 



(r2) = / 

Jv 



= I dvs 

V 



(G23) |'«(r3)a„(r3) (A3) 



In both cases, Eq. (A2| and (A3), the incoming inten- 



sity Cin may be of type Ci or Cc (coherent mode) . (In 
the latter case, r and k must be replaced by Tc and Kc, 
respectively.) Ther efore, th e transport equation for Ci is 
obtained from Eqs. ( |A2|A3| by setting Ci = C(^)-hC(°+'=) 
and inserting Gin — Ci + Cc- Similar considerations 
can be repeated for the calculation of the coherent mode 
Cc- The resulting transport equations, graphically rep- 



resented in Fig. 11 read as follows 
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Cc{v) 
Ci(r) 



(A4) 



dr'|(G(r,r'))r 



k{v')C,{v') + k,{v')C,{v')+ 



+ (r(r')Ci(r') + T,(r')Ce(r') 



(A5) 



|r'-r|_ 



Furthermore, 



with q = Ul + kD and ri = r' + pi 
Tz' and z„ are defined as in Eqs. (26 1 and (|28|. Note that 



in our notation for the incoming or outgoing waves, we 
have suppressed the position dependence of the refractive 
index. The latter must be averaged over the correspond- 
ing path, e.g. exp(z/cnc-z) '■— exp[ifc dz' nc{v z')]- 



3. Equations for C2 

Concerning the crossed intensity of type C2 , the equa- 
tion turns out to be more comphcated. For example, 
it may happen that two nonlinear propagation events. 
Fig. [6]je) and (f), occur between two scattering events, 
and, furthermore, all the cases must be distinguished 



where the incoming solid or dashed arrows do or do not 
originate from the coherent mode. Graphically, the result 
is shown in Fig. [T2j 

In order to shorten the notation, we split the equa- 
tion into three parts. Here, Co ("o" for "outgoing") en- 
compasses all diagrams, where the dashed arrow origi- 
nates from the incoming laser and, correspondingly, the 
solid arrow points towards the detector. This quantity 
finally determines the coherent backscattering cone F,^, 
see Eq. (A9). The quantity Cp ("p" for "propagating") 
contains the part of C2, which originates from the non- 
linear propagation event Fig. [6]jf ) occurring between two 
scattering events. 

The mathematical expression of Fig. [12] yields the cor- 
responding transport equations as follows: 



C2(r) 



k{v)C^{v) + k,{v)C,{v) + T*{v)Cp{v) + K*{v) [ rfr'|(G(r,r'))|'C2(r') 

Jv 

C2{v') + k{v')C^{v') + Ko{v')C,{v')+ 



C,(r) = / dv'\{G{v,v'))t 



-(r(r')Ci(r') + re(r')Ce(r')) / dp.e'P-" K{v^) 



Co{v) 



ik{n~-n*)z 



^«(r) + ^:(r)) C\{v) + r;(r)Cp(r) + <(r) 



dp2e-P'/'K{r2) 
dr'\{G{r,r'))fc,{r') 



with n = r' + pi ^r^^ and ra 



anu 1-2 =r-Hp2|7=^7|. 

Finally, the height of the coherent backscattering cone 
reads: 



(A6) 



(A7) 



(A8) 



rc(kD) = / 

Jv 



dr 



4. Simplified version 

It remains to be shown how the simplified version, see 
Eqs. ( |45p8l ), is derived from Eqs. ( [A4]|X9| . For this pur- 

rtnsp wp flssnmp fVij^f, wVipn fViP rinsifinn in 



pose, we assume that, when varying the position inside 
"■^(r) (A9) 

the sample volume on a length scale comparable to the 
extinction path length, the average intensity remains ap- 
proximately constant. This is the case for a medium 



19 



K 



C2(r)= C2: 



+ 



a, b 



2<., 



+ 



1 1 



K(y) 

A 



Figure 12: Graphical representation of the transport equa- 
tions for the crossed intensity of type C2, Eqs. ( A6|A8 1. 



with large optical thickness, except for the regions close 
to the boundary. Using this assumption, and setting 
^(max) _ ^ (which is again valid far away from the 
boundary) , we simplify the integral over pi in Eq. ( A5 ) 
as follows: 



dpiexp(-pi/£)X(ri) ~ if(r')^(r') 



(AlO) 



This yields the simplified equation for Ci, see Eq. (46 1 



with the cooperon cross section a, Eq. (43 1. In the equa- 



tion for the coherent mode Cc, Eq. (A4|, we neglect the 
term containing t^, and obtain Eq. (451. This is justi- 



fied for a relatively weak nonlinearity, where a nonlinear 
event is unlikely to occur on the length scale £ defined by 
linear scattering. 

Using the same argument, we neglect the term con- 
taining T* at the end of Eq. (A6). Thereby, the Co-part 



drops out of the equation for C2, whereas the Cp-part is 
approximated by 

C,(r) J^dr'\{G{r,r'))f[c2{r') + a{r')C^{r') 

+a,{v')C,{v'))K{r)£{v)^ (All) 
if(r)£(r)[Ci(r) + C'2(r)] (A12) 

where Eq. pG] ) was used, and we defined 

C2(r)= / dr'|(G(r,r'))|'C2(r') (A13) 
Jv 



Finally, t^ 



Now, transforming Eq. ( A6 ) into an equation for C2 
yields Eq. (|47l 



;he bistatic coefficient, Eq. (48 1 
is obtained from Eqs. ( A8|A9 1 under neglect of the last 
term (the one containing Tc) in Eq. (AS I. 

In summary, the simplified transport equations. 



Eqs. ( 45][48 1, are valid in the case of large optical thick- 
ness (6 ^ 1) and weak nonlinearity (such that nonlinear 
events close to the boundary of the scattering volume can 
be neglected) . Note that these two conditions are consis- 
tent with one another in the following sense: for larger 
optical thickness, a weaker nonlinearity is required to in- 
duce a comparable effect on the coherent backscattering 
cone. 



Appendix B: PROOF OF THE RELATION ([SSj) 
BETWEEN COOPERON AND DIFFUSON 
CROSS SECTIONS 



In this appendix, we will proof the equations 

Id d 



Re k\ 



/c=0 



'Re{T*K}L, 



2 did 



la=0 
d 1 



2^ did £ 



(Bl) 
(B2) 



from which Eq. ( 53 1 results according to the definition of 
(T, see Eq. (Ilif. 



To show Eq. (Bl I, we start with Eq. (|38| for k, take its 
real part, and express the expectation value (...) as an 
integral over the intensity distribution function P{I), see 
Eq. (17 1. Finally, we insert / = Eg{I), and express the 
field derivatives d/dE and d/dE* as intensity derivatives 
d/dl. Thereby, we obtain: 

f°° d d 

Rek = N dIP{I)-I^-\g{I)\^ (B3) 

= dl {p'{I) + ^-P"{I)^ I\g{I)\' (B4) 



where we have used partial integration. Since Eq. (Bll 
is evaluated at 1^ = 0, we take P{I) = cxp{— I / Id) /Id, 
see Eq. ^ with = 0. 



In a similar way, we obtain 

(7W = -AA J dIP\I)I\g{I)\' 



From this, Eq. (Bl I follows with help of 



|Ap'(z)^P'(/) + ^P"(/) 



(B5) 



(B6) 



valid for P(/) = exp(-///<j)//rf. 

To show Eq. (|B2), we first note K ^ Af j dIP{I)I\g\'^ , 
see Eq. and~P(/) = -IdP'{I). Using Eq. (|B5|, we 

see that Eq. (B2| is equivalent to 



Rer* 



1 d 1 

2~drd~£ 



(B7) 



/e=0 
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Now, we rewrite Eq. (41 1 as 



2fc 



(B8) 



dl {IP"'{I) + 2P"(/)) Ig{I) (B9) 



2fc Jo 

On the other hand, Eqs. ( |20|2ip yield 

dIP'{I)Iluig{I) (Bll) 



Finally, Eq. (|B7| results from Eqs. ( |B9|B11[ ) and 

d 



dl, 



P'{I) = IP"'{I) + 2P"(/) 



(B12) 



for P{I) = exp(-///d)//d. 



We note that Eq. (Bll refers to nonlinear scattering, 
and Eq. (B2) to nonlinear propagation diagrams. Since 
the equality holds independently in both cases, the above 
derivation, explicitly performed for the case of nonlinear 
point scatterers, see Eq. ([2]), can be generalized to other 
nonlinear models, where scattering and propagation have 
different physical origin. In particular, this is the case in 
our second model, i.e. linear scatterers embedded in a 
homogeneous Kerr nonlinearity, see Eq. (54 1. 
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